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Abstract 

We derive the conditions that the coupling constants of the Generalized Thirring 
Model have to satisfy in order for the model to admit an infinite number of commuting 
classical conserved quantities. Our treatment uses the bosonized version of the model, with 
periodic boundary conditions imposed on the space coordinate. Some explicit examples 
that satisfy these conditions are discussed. We show that, with a different set of boundary 
conditions, there exist additional conserved quantities, and we find the Poisson Bracket 
algebra satisfied by them. 
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1. Introduction 

The purpose of this paper is investigate the conditions under which the Generalized 
Thirring Model becomes integrable. The Generalized Thirring Model is the field theory 
of massless fermions in two dimension, interacting through the most general four fermi 
interaction compatible with Lorentz invariance [1]] . By integrability, we mean the existence 
of an infinite number of conserved and commuting dynamical variables. In this paper, we 
will not address the question whether these are sufficient in number to make the model 
truly solvable. The bosonized form of the Thirring model will be our starting point, and 
our treatment from that point on will be purely classical. It should, however, be noted that 
this is better than treating the original fermionic model classically, since bosonization does 
capture some of the quantum nature of the model. From the bosonized Lagrangian, we wish 
to extract a Lax pair depending on a spectral parameter. For this purpose, we first write the 
equations of motion in a suggestive form as flatness conditions for two vector fields, and we 
demand the existence of another flat vector field that interpolates between these two. This 
results in equations involving the coupling constants which we call the first integrability 
condition. In general, these equations are overdetermined and they do not have solutions 
depending on a continuous parameter. In section 3, we discuss four exceptional cases when 
such a solution exists. The first case is the model with maximum internal symmetry and 
a single coupling constant. The second case is a simple generalization of the first one to 
a product group with two different coupling constants. The third example is an 577(2) 
model broken down to U(l). The last example is a model based on symmetric spaces. For 
each of these examples, there exists a Lax pair depending on a spectral parameter. 

Given the Lax pair, in a standard fashion, one can construct conserved quantities in 
terms of a path ordered product. It is, however, necessary to specify boundary conditions. 
The simplest boundary condition is the periodic one, with the space coordinate compacti- 
fied into a circle. Taking the trace of the path ordered product and expanding in powers of 
the spectral parameter yields an infinite number of conserved "charges" . In section 4, we 
compute the Poisson brackets of these charges and derive the conditions so that it vanishes. 
This is then our second integrability condition. In the four examples discussed earlier, the 
second integrability condition is automatically satisfied, although we are unable to prove 
in general that the second condition follows from the first. Another natural boundary 
condition is the open one: The space interval is from — oo to +oo and the fields vanish 
at ±oo. In this case, additional integrals of motion can be constructed by considering the 
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matrix elements of the path ordered product, instead of just the trace. In section 5, we 
compute the Poisson brackets of these additional integrals of motion and find that they 
satisfy a non-linear algebra. We also compare this algebra to a simpler algebra derived in 
a somewhat similar case of the principal chiral model . The last section summarizes our 
conclusions. 



2. Lax Pair Formulation 

We begin this section by recalling the definition of the model, the bosonized version 
of it and its equations of motion and symmetries (see ]1J and for detailed analysis). 
Lightcone variables will be used throughout: x_|_ will serve as time and X- as space. 

The parity violating generalized Thirring model is given by the action 

S = J d 2 x{m^d^ - {G~ x ) ab y R t a V R *l*6*l), (2.1) 

where R and L refer to the right and left chiral components of and t a are the generators 
of the Lie algebra Q in some representation: 

[ta, tb] = if ab C t c . 

We reserve the notation t a for the adjoint representation. The coupling constant (G~ 1 ) ab 
is an invertible not necessarily symmetric matrix, resulting in parity violation. In one 
version of bosonization Q [|J , this gives 

So = W(g) + W{h~ l ) [ d 2 xG ab (ig- 1 d + g) a (ih- 1 d_h) b , (2.2) 

Zn J 

where X a stands for Tr(t a X), g and h are group elements expressed in the adjoint repre- 
sentation, and 

1 7T ~ 

G a b = Ti^ab — -^-Gab, (2.3) 

where K a b is the Cartan-Killing metric n a b = Tr(t a tjj), n is the number of fermion flavors 
shifted by the dual Coxeter number and W is the WZW action 

W(g) = ^ (/ d'xTrid^g-^g) + \ j Tr ({g^dgf)^ . (2.4) 
The equations of motion are equivalent to conservation of two currents: 

d + J_ = d_J + = 0, (2.5) 
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where 

J + = i— (-d+hh- 1 + ht a h- 1 (2G T ) a b (g- 1 d + g) b ) , 

(2 6) 

J_ = z£- (-d-gg- 1 + gt a g-\2G) a b {h- l d-h) b ) . 

Next we notice that the model is invariant under the transformation h — > u+(x+)h, with 
J-f transforming like a gauge field. Using this transformation, J+, which only depends on 
x+, (see (|2.5|) ), can be set equal to zero. This gives us a special solution to the equations of 
motion; the general solution is obtained by applying the inverse transformationB. For this 
special solution, treating x+ as time, we will search for an infinite set of conserved quanti- 
ties, which assure the integrability of the model. The general solution is also integrable by 
virtue of the transformation introduced above. The situation is similar to what happens 
in the WZW model |3j . The conserved quantities are derived for this special solution with 
J_i_ = 0. Define now 

VI = (ih~ 1 d ± h) a , (2.7) 
W$ = (2G T )- la b V b , = (2G) a b V*. (2.8) 
In terms of these variables, the equations of motion now read 

d+V- - d-V+ - i[V+ , V- ] = 0, 

(2.9) 

d+W_ - d_W + - i[W+, W-] = 0. 

These flatness conditions for the vector fields V and W are similar to the zero curvature 
condition of integrable systems. What is missing is a spectral parameter dependence that 
will provide by power expansion an infinite number of conserved currents. We will now 
derive conditions for a zero curvature with a spectral parameter dependence (a Lax pair) 
to exist. This will be the first integrability condition. The idea is to find a one parameter 
family of matrices that interpolate between the equations of motion. It is convenient for 
later use to work with a special linear combination of the V± connection: 

Ml = {Hl{2G T )- 1 ) a b V b , H+ = l- 4GG T , 
M a _ = (Hl) a b V b , iZ_ = 1 - 4G T G. 

Define now an interpolating connection as follows: 

Bl(x;\) = N£ b (\)Mi(x), (2.10) 



1 We thank Bogdan Morariu for clarifying this point for us. 



3 



where N± b are constants, to be determined as functions of the spectral parameter A. The 
vector field B must satisfy the zero curvature condition 

d+B_ -d_B + -i[B+,B_] = 0, (2.11) 

with boundary conditions 

SJ(x;A = A ) = K t a (x), 

(2.12) 

B^(x]X = Ai) = Wl(x). 
In order to proceed we rewrite, after some manipulation, the equations of motion as follows: 

d-Ml + A a _MlM q _ = 0, 

Pq ' (2.13) 
d+M* + A% vq M v _M\ = 0, 

where 

A ±abc = -(^/(tf-^^G^)^ (2.14) 

and for convenience we defined G+ = 2G T and G- = 2G. Note that A± a i, c is antisymmetric 
in its first two indices. The zero curvature reads 

N* h d + M b _ - Nl h d_M h + + f a bc N b +v N^ q MlM q _ = 0. (2.15) 

Using equations (|2.13|) and equating the coefficients of M+Ml to zero we get the first 
integrability condition 

N a +h A b _ pq - N a _ h A b +qp + f a bc N b +p N c _ q = (2.16) 

This condition gives us a Lax pair and an infinite number of resulting conservation laws 
(see the section 4). One must also show that these conserved quantities are mutually 
commuting; that is, their Poisson brackets vanish. This will be the second integrability 
condition, derived in section 4. 

Equations (|2.16|) are in general an overdetermined algebraic system with (dimG) 3 
equations for 2(dimG) 2 variables N± b . In fact, since these equations are nonlinear, this 
counting is misleading. What we have actually is a system of polynomial equations for 
the variables N± h . The locus of the polynomials defines an algebraic variety M. and the 
condition of integrability is not that Ai ^ {0} (this is guaranteed since we always have 
two solutions when B is equal V or W) but that dim.M > 1 in order to have a spectral 
parameter. 

Although the system is overdetermined and there are no parametric solutions for 
a generic coupling constants G a b, there are special interactions for which the model is 
integrable. This is the subject of the next section. 
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3. Solutions of the First Integrability Condition 

In this section we will construct some solutions to the integrability condition, proving 
thus, the existence in those models, of an infinite number of conserved currents. In all 
examples we will make use of the diagonal ansatz. We assume that G and N± are diagonal 
matrices: 

2G\ = g b 5% 



Then A± abc = A abc , 



N a ± b = KS% 



a 9a9b 9c n / q -. \ 

A-abc — 7~ i „ i 'Z~TJabct 



and the integrability condition (|2.16|) reduces to 

n+A abc - n~A acb + f a b c n£n~ = 0, 

or, using (Q , 

/ -f 9a9b ~ 9c _ 9a9c ~ 9b + _\ ~ _ _ 

(1 - 91) 2 (1 - b ) 2 (1 - 0c) 2 (1 - 91) 2 (1 - 9b) 2 (1 - 9c) 2 



(3.2) 



Note that there is no implied sum here. 

Example 1: 2G = gl. (The Symmetric Case) 

In this case, N± = n^l, and we have one equation for two variables, 



2 2 

n + ^ + n~ ^-+n + n~=0, (3.3) 

(1-^)1 (l-g*)i 



with a one parameter family of solutions. 
Example 2: 2G = g±l <g> (72I 

This just gives two copies of the above equation ( |3.3|) . 
Example 3: SU(2) with £7(1) symmetry. 2G a b = g b S a b with #i = #2 7^ 93- 
This suggests the ansatz, iV| b = n^d^ with = ^ nf. Then there are three 
equations for four variables, 

n t(gi - gs) + nligm - 91) + (1 - g\)0- - gi^nfn^ = 0, 

"3 (0301 - gi) + ^3 toi - 0i) + (1 - 0?)(1 - gl)^n+n^ = 0, (3.4) 
nt{9i93 ~ 0i) + %(0i - 03) + (1 - 0i)(l - 03)^3"% = °- 
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The one parameter family of solutions is given by, 



t = {mTT^W (2{93 ~ 9?) + ^ 9l(1 " 93) + 91(1 " 93,A) ) 

1 

l" = ( (1 + gs)(l-g^)2 ^ 93 ~ 9 ^ + \ 91 ^ ~ ^ + ~ ^ 3)A )) 



n 



l 3 = 7i ^(tes-sD + T^ 1 -^)) 

3 = 7j iTTTi -^{(93- gl) + Xgi(l - g 3 )) ■. 



(3.5) 



(l-^)ni-^) 



where A = ^qr is the free parameter with range ^- < A < gi. Note that the case where we 
also have (72 7^ #1 gives six equations for six variables. 
Example 4: Symmetric spaces 

Let F be a simple group with a subgroup H. Then the Lie algebra T can be decom- 
posed into the Lie algebra TC and its orthogonal complement /C, which generates the coset 
F/H. This coset space is a symmetric space if [/C,/C] C T~i. In what follows we will label 
the generators of H with greek indices (i.e. A a ), and the generators of F/H with latin 
indices (i.e. A a ), and when we don't want to specify between them, we'll use dotted latin 
indices (i.e. A a ). 

To have a coset (symmetric) space in the present model we choose 

2 <n = 9iS%, 

with g a = g and g a = 1. This assures that the currents in 7i are set to zero, resulting in 
a coset model. We cannot use ( [2.16|) because A± a t, c is not defined in this case (H± is not 



invertible), but instead, we have to use (|2.9|) and (|2.11|) 



d+V* - d_V? + f\.VlVi = 0, 



+ ^ J bc v + 

a 
be 

a 
b 

The choice of G leads to the ansatz: 



d+Wl - d-Wi + f\.WlWi = 0, (3.6) 
d + B h _ -d_Bi + f\.B h + B d _ = 0. 



B a ± (x;\)=r±(\)Vi(x), 
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with r^(X) = r ± (A) and 7™(A) = 1. Then, we can rewrite (|3.6| ) as 



(i - gl)d + v. a + f, b .(i - 9 -f^)v +b v. c = o, 

(1 - \)d-V + a + /^(l " -7)^ = °' (3 - 7) 
yd ya% 

r-d+V-a ~ r+d.V +a + f abc r+r7 V+b V. c = 0. 
Note that a is fixed (no implied sum). Now if d is in Tl, then the first two equations in 
( p.6|) are the same, and from the third one we get the condition r + (A)r~(A) = 1. The 
fact that the first two equations become the same is a sign of gauge invariance. To fix a 
gauge we choose V- a = 0. Then in this gauge, if a is in /C, the first two equations in Q3.7Q 
solve the third one without further conditions. Hence, for symmetric spaces, we get one 
equation for two variables 

r+(A)r"(A) = 1. (3.8) 

4. The Poisson Bracket — Periodic Boundary Condition 

Having a Lax pair at hand, we can construct conserved quantities in the time variable 
x+, if we also impose periodicity in the space coordinate x_. We first define the following 
quantity: 

U(x +I x-;X) =Pe Jx - ~ (4.1) 

and take B± periodic in X- with period 2-k. The integral goes between X- and X- + 2n ~ 
X-. This quantity satisfies the equation 

d.U -i[B_,U] = 0, (4.2) 

which can be taken as the definition @ of path ordering in (|4.1| ). More importantly, the 
trace of this matrix, U(\) = Tr U(x+,x-; A), is conserved: 

d+U = i ( Tt(d+BM)dx'_ 

J X- 



' X - 
f X- 



(Bianchi id.) = i ! Tr Ud-B+ + i[B+, U)dx'_ 

(4.3) 

(integration by parts) = — % \ Tr (B+(d-U — i[B-,U]))dx'_ 

J X 



(by O) = 0. 



2 This is the same as denning the path ordered exponential with the "upper" limit fixed, i.e., 

f b dA 

Pe-ia = 1 + J a dxA(x) + J a dx f x dx A{x)A{x) + • • •. Note that this is the opposite of the 



usual definition where the "lower" limit is fixed. 
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We see that U is a conserved quantity and upon expanding it in powers of A we get infinite 
number of conserved currents. Another possible way of of increasing the list of conserved 
quantities is to choose B to live in different representations of the Lie algebra. Since U 
does not depend on x+, we drop its dependence from U as well as the subindex — from 
X-. It is convenient to think on the + direction as "time" and the — direction as "space". 
The next step is to find the algebra of the conserved currents. In particular if they all 
commute then we have infinite number of conserved quantities in involution, which is the 
trademark of integrability. We want to calculate the Poisson brackets of U(X) with U(/j,): 

MA) , uv) } = J d**» ™_$ x) jffe <*--fr 

= - j dxdyU a (x;X)U b (y-fi){B a _(x;X),B b _(y,^} (4-4) 
dxdyU a {x;\)U b {yiriN° e {\)N± d (ri{Ml{x),Mi{y)y 



The Poisson brackets for the M's where calculated in fill : 



47T 47T 

{M c _(x), M d (y)} = K cd 8'(x -y) + —Fl d k MHx)8(x - y) 

n n 

+ —e(x-y)E cd kl (x,y)M^x)M l _(y), 
n 



where 



(4.5) 



E c U^y) = A^A^(Pei: A - M -^ dx ') rs , (A. c ) ab = A_ abc , 

F +abc = (HZ l2 )a r (HZ h ) b s (HZ l2 )cfrst - SiHZ^G^JiHZ^G^iHZ^G^^f, 



(4.6) 



Notice that the Poisson brackets include a non local term! It is clear that the only way to 
get a local algebra for the Poisson brackets for our conserved fields is that the coefficient 
of e(x — y) is a derivative such that by integration by parts the derivative acts on e to give 
a 8 function. For a generic coupling this can not happen, but exactly for the models that 



2tt 
n 



satisfy the integrability condition this is true! 

{U(X) , U{jl)} non-local = 

= *L f dx^ a (x;A)^(y;^)Ar« c (A)^ d (^) e (x-y)i? c i z (x,y)M^(x)Mi(y) 
n J 

= ^1 J dxU a {x;X)N1 c {X)A* k r M*{x) 

J dy e (x-y)^(y^)(Ar^(^)^ s )(Pe/>- feM - (x,)dl ') rs Mi^) (4-7) 
dxU a (x; A)r c (A)A + VM^x) 

dye(x - vWttoMNXaMA^ + f 6 tu iV|>)iV^(/i)) 

x {pJ>- kM - {x ' )dx ') r M l _(y), 

where the last step follows from the integrability condition, ( |2.16| ) . Using the following 
identities 

d y Ut(y^) = -f b tu B u _(y^)U b (y^) = -f b tu N^)M l _(y)U b (y^), 

we can write now the integrand of the y integral as a total derivative 

Ut(v,») {NlMA^+f^^NlM) {PeJI A -- M ^ x "> dx ') r M'_(y) = 

-a,{u b { v ^)NXM(pJ>-^)"- (49) 



The contribution of the non-local piece is therefore 

{[/(A) , U(v)}non-iocai = ^- J dxU a (x; A) U b (x; fi) AT« c ( A) N b _ d (ji) A'fMl (x) 

= 2 1 j dxU a (x;\)U b (x;ri (N^ c (X)N^ d (fi)A_^ l d - N^ c (X)N^_ d (fj)A d l c ) M l _(x), 

(4.10) 

where in the second step we (anti) symmetrized the right hand side so that the change in 
sign under A <-> \i is manifest. The contributions from the local terms are easily evaluated: 

{17(A) , U(jj)} s = ~J dxUa&Xpt&riN^NlMFfiMlix), (4.11) 



and 



27T f 

{U(X) , U(fjL)}s> = — / dxU a (x; X)U b (x; n) 

n J (4.12) 



Putting everything together we get 

{17(A) , U(n)} = J dxU a (x; X)U b (x; fj)J ab i(X, n)M l _(x), (4.13) 



where 

J a \(X, y) =N a _ c {X)N h +d ^)A< l d - Nl c (X)N b _MA^ c - 2N- c (X)N b _ d ( fJ ,)F^ d 



(4.14) 



The right hand side of equation fl4.13| ) is zero, if the integrand is a total derivative. This 
motivates us to try the ansatz 

d x (U a (x;X)U b (x;v)C ab (X,Li)) = U a (x; X)U b (x; ^J a \(X, ^)M l _(x). (4.15) 



The left hand side of this equation can be computed using the identities of ( |4.8| ) . This 
result in the second integrability condition 

J o6 Am) = r cd N^(X)C cb (X^) + f b cd N^)C ac (X^). (4.16) 

In the diagonal ansatz, C has the form 

C a6 (A, fi) = C a (A,//K 6 . 

For models for which a C ab that satisfies this equation can be found, {U (A), U (fx)} = 0, and 
the conserved quantities are in involution. We will now check the models that satisfy the 
first integrability condition against this second integrability condition. In all examples we 
have, the result is the same: The second integrability condition is automatically satisfied. 
However, we do not know whether in general the second condition follows from the first 
one. 

In example 1, the symmetric model, we take C a (A, fi) = C(A, fj,), and then we have one 
equation for one unknown, C(X,fi). The condition ( 14.160 is satisfied, and the conserved 
quantities are commutative. The generalization to example 2 is trivial. In case of example 
3, SU(2) with C ab (X,fi) = C a (X,y)K ab and with C X {X^) = C 2 (A 5/ u) ^ C 3 (A,^), initially 
we have three equations for two unknowns, but using the fact that C a (A, /u) = —C a (p, A) 
we see that there are only two equations. So, in this case also, equation ( f4.16| ) can be 
solved, and the model is integrable. We have not checked example 4 in detail, but we 
suspect that the result is the same. 
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5. The Poisson Bracket — Open Boundary Condition 

In this section we will show that the structure of the model is richer with open bound- 
ary conditions. We consider the open interval from — oo to +00, with fields vanishing at 
±00. With these conditions, it is possible to construct more integrals of motion than with 
the periodic one. These integrals do not commute in general and it is the aim of this 
section to calculate the algebra they generate. A subset of these integrals are the integrals 
we saw in the periodic case, and they commute, insuring, thus, the integrability of the 
model. 

The calculation of the Poisson brackets is the same as the periodic case, the only 
difference will come from the boundary. Defining 

U x {x,y)=Pe- l S>- {x '' X)dx \ (5.1) 

then, any matrix element U ab (X) = U% b (— 00, 00), and not just the trace, is conserved, 
provided that S + (±oo, A) = 0. This follows from manipulations similar to Ql.3|) , noticing 
that the boundary terms coming from ±00 can be dropped. Of course, B- also must 
vanish sufficiently fast at ±00 for the integral to exist. For the sake of simplicity, we will 
also take B in the adjoint representation. 
Define now 

U ab c (x; A) = (W A (-oo, x)tJA x {x, oo)) ab , (5.2) 

then 

/XTTaa' ( \ \ STI bb '(ll) 
^ ^ ( ; ; A)^ ( ^) { ^^ A) ^- (y;/A)} 

= - J dxdyU aa c(x;\)U b %^)Nl k (\)N^){M k _(x),^ 

(5.3) 

Being careful to take into account the boundary contribution from the non-local term in 
the Poisson brackets of the M's, we get 

{U aa '(X) , U bb 'U)}=— [ dxU a i(x;\)U bb ' d (x;^)J c i(\^)M e _{x) 

n J 

+ ^(f/(A)g aa '(^(^)) 6b V^(0)iV^(A)iV^(^) (5.4) 
-^(^(A)) aa '([/(^)t d ) 6b V fcZ (0)iV^(A)iV^^), 

11 



where J ab c is defined in ( |4.14|) , and U ab (0) = U ab (X = A ). If J ab c also satisfies the second 
integrability condition ( f4.16| ) then the integrand can be written as a derivative. The only 
contributions come from the boundary and the final result is 



{u aa' {x) j ?7 66 '(/i)} = —((U(\)t c ) aa '(U(fl)td) bb C Cd (\,fJ,) 

-(t c U(X)) aa ' \t d U{y)) W 'C cd {\^) 



(5-5) 

+ {U(X)t c ) aa {t d U{y)) bb U lk (0)N c -k(mUri 
-(t c U(\)) aa \u^)t d ) bb 'u k \0)Nl k (X)N^)). 



This algebra, in the symmetric case, can be compared with the algebra obtained without 
bosonizing the fermions @ . Both the appearance of cubic terms and the more complicated 
spectral dependence of the "classical" r-matrix C a6 (A,/i) are quantum effects due to the 
process of bosonization. 

Equation ( |5.5[ ) exhibits a closed algebra for the matrix elements U ab (X); iV's and C's 
are the corresponding "structure constants" that depend on the particular model under 
consideration. The algebra is clearly non-linear, with quadratic and cubic terms in U 
appearing on the right hand side. We have not succeeded in identifying it with any well- 
known algebra, although it may possibly bear some relation to the W algebras ||. It is 
also the generalization of the algebra Liischer and Pohlmeyer || found for the special case 
of the fundamental representation of SU(2). Finally we would like to comment on the 
relation of this algebra to the affine algebra of currents found in the principal chiral model 
0. In the symmetric case (Example 1 of section 3), the Thirring model bears a great 



resemblance to the principal chiral model; for example, the equations of motion (|2.9|) in 
this case reduce to 

gd+V- +d_V + = 0, 

(5.6) 

d+V- - d-V+ - i[V+, V-] = 0. 

For g = 1, these are identical to the equations of motion of the principal chiral model. Not 
surprisingly, the Lax pair and the conserved quantities are in one to one correspondence, 
with obvious modifications in case g ^ 1. However, the algebra (|5.5|) is quite different 
from the affine algebra found in @. This is because in ||, the Poisson brackets of the 
conserved quantities were derived from the transformations they generate on the field 
variables, whereas we have used the standard Poisson structure given by the Lagrangian. 



The two Poisson structures differ in this well as in the principal chiral model |TIJ . 
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6. Conclusion 

In this paper we derived very general conditions (equations ( |2.16|) and (|4.16|) and 
Q4.14T) ) that the coupling constant of the Generalized Thirring Model should satisfy to 
be integrable, in the sense of having an infinite number of conserved quantities. Some 
solutions were found for special cases. The algebra generated by the conserved currents 
was calculated for both periodic and open boundary conditions. In the periodic case the 
currents are in involution while in the open case we found a non linear algebra which has 
an Abelian subalgebra of conserved currents in involution. 

The directions for further research are numerous. One may, for example, add fermions 
to the bosonized theory and study supersymmetric models. Or, one can try to generalize 



to a space dependent metric |1JJ . There is of course the problem of finding new solutions 
to the overdetermined algebraic systems eqs. ( |2.1(3| ), Qll6| ) and ( [4.14| ). It is important 
to understand the non linear algebra between the conserved quantities that emerges in 
the open boundary case.. We would also like to quantize these models and to study the 
quantum algebraic structure. As a first check we would like to see if the relation between 
the coupling constants that ensures integrability holds along the renormalization group 
flow. For the solutions we have studied in this paper, this result follows rather trivially 
from symmetry considerations. Finally there are the questions of the relation to irrational 
conformal field theories (see |12| and references therein), and to integrable perturbations 
of conformal field theories. We hope to address some of these questions in the future. 
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